Lost and found: The missing diabolical points in the Fe 8 molecular magnet 
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Certain diabolical points in the tunneling spectrum of the single-molecule magnet Feg were pre- 
viously believed to be have been eliminated as a result of a weak fourth-order anisotropy. As shown 
by Bruno, this is not so, and the points are only displaced in the magnetic field space along the 
medium anisotropy direction. The previously missing points are numerically located by following 
the lines of the Berry curvature. The importance of an experimental search for these rediscovered 
points is discussed. 

PACS numbers: 75.50.Xx, 03.65.Vf, 03.65.Db 



The purpose of this note is to report on a numerical 
search of certain diabolical points (DP's) in the energy 
spectrum of the single-molecule magnet (also known as a 
molecular magnet or nanomagnet) Fes that were earlier 
believed to be missing, but are in fact not so Sev- 
eral other DP's have been seen experimentally in Feg Q, 
and their observation provides the best evidence of spin 
orientation tunneling between deep levels in all single- 
molecule magnets studied to date. Observation of even 
some of the missing DP's would strengthen our under- 
standing of Fes substantially. 

For a system whose Hamiltonian depends on some pa- 
rameter, a DP is a point in parameter space where two 
(or more) energy levels are degenerate [3j. In Feg, the 
parameter is the static applied magnetic field, and the 
locations of the DP's so far observed (as well as many 
other experimental measurements) are well described by 
the anisotropy Hamiltonian, 



H = k x Jl + k 2 jl - C{J% + Ji) - gn B 3 ■ H. 



(1) 



Here J = (J x , J y , J z ) is the spin, H is an external mag- 
netic field, g is a g- factor, fis is the Bohr magneton, and 
ki, &2, and C are anisotropy coefficients. Experimentally, 
J = 10, h ~ 0.338 K, k 2 ^ 0.246 K, C ~ 29 /jK, and 
g ~ 2. The DP's can be understood as arising when 
tunneling between two states with (nearly) oppositely 
oriented magnetic moment is quenched because of de- 
structive interference between instantons (spin tunneling 
trajectories) [3, HJ. 

The model Hamiltonian (fT]) was first analyzed in 
1993 [4] with C = 0, and it was found that for ground 
state tunneling it had 10 DP's along the positive H x axis, 
corresponding to J = 10. In reality only 4 DP's are 
seen [2j, which was explained in Ref. [6J as follows. When 
C ^ 0, we get two new (but noninterfering) instantons, 
which are discontinuous at the end points. One of these 
instantons has the least action when H x exceeds a cer- 
tain value H* , and since this instanton has no interfering 
partner, there are no more DP's for H x > H* . For the 
values of the anisotropy coefficients quoted above, H* lies 
just beyond the location of the fourth DP, which explains 
why the last six DP's are not seen in experiment or by 



direct numerical diagnolization of Eq. (TTJ). We can fur- 
ther check this picture by decreasing C so as to make H* 
larger. It is seen that as this happens, one successively 
sees six and then eight, and finally all ten DP's, all in 
accord with direct numerical diagonalization results @ . 

However, as shown by Bruno 1], the above picture, 
though correct, is incomplete. For any energy level, the 
sum of the Chern numbers for all DPs involving that 
level is a topological invariant as parameters like ki, ki, 
or C are varied. Since DP's in any system are gener- 
ically simple, we expect this to be so in Fes also, and 
the Chern number for any one DP should be ±1 whether 
C = or C 7^ 0. Hence the six missing DP's must 
be present elsewhere in magnetic field space. A similar 
conclusion applies to the DP's associated with tunneling 
between other pairs of levels 0] ■ For tunneling between 
the ground states, the DP's merely move off the x axis 
into the xy plane. For the higher energy levels, they move 
off the xz plane into the full three dimensional H space. 
This point can also be understood by noting that for a 
system with purely four- fold symmetry (ki = hi, C ^ 0), 
the ground state DP's lie on the ±x ± y axes, while for 
the excited states they lie in the planes formed by these 
axes and the z axis [|| . When both two- fold and four- fold 
anisotropies are present (ki ^ ki, C ^ 0), it is then not 
surprising that the location of some of the DP's is also 
intermediate [l|. 

Observation of these rediscovered DPs would be in- 
teresting in itself, and also provide an important test of 
the validity of the model (JTj) vs. other models that 
add extra 6th and 8th order anisotropies because the lo- 
cation of the DPs is very sensitive to the higher order 
anisotropies. With this motivation, we have undertaken 
a search for the DP's for the ground state and some of the 
excited states. We stress that the key insight that these 
points should exist in the first place is due to Bruno, and 
our contribution is only to find their specific locations. 
Neverthless, finding them is not without challenge as we 
discuss next. 

A direct search for the DPs by numerical minimiza- 
tion of the energy differences fails because the energy 
surface is like a golf course with rolling hills on which 
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TABLE I: Locations of selected diabolical points 
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FIG. 1: Diabolical points of Fes corresponding to tunneling 
from the Zeeman level m = — 10 to m — 10 (layer 1), to m — 9 
(layer 2), and to m — 8 (layer 3), all projected onto the xy 
plane. H c = 2k\J / g^B- 



the DPs are the holes. Because the holes are so local- 
ized, unless one starts close to one of them by luck, any 
numerical algorithm will in general simply head for the 
valleys of the course and miss the holes entirely. Be- 
cause H is not real for general H, we also cannot use 
the method of Ref. |8j, which is to corral the DPs by 
using the Herzberg and Longuet-Higgins theorem |9j to 
find and successively bisect a sign- reversing circuit. Wc 
therefore proceed as follows. Let us denote the eigen- 
states and eigenvalues of Eq. ([T]) for fixed H by |n(H)) 
and E n (H), n = 1, 2,..., 21, and order them so that 
E n > E n -i for every H. Except at degeneracies (the 
DPs), the Berry curvature for the nth level is defined by 
111 



B, 



j v-* (n\V H H\n') x (n'\V H H\n) 
Im 1^ 71? — • ( 2 ) 



(E n , - E n y 



Further, the Chern number associated with a degeneracy 
is given by 



Q, 
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B„ • d s, 



(3) 



where the integral is independent of the choice of the 
surface S, as long as it encloses the degeneracy, since 
Vh • B n = away from the degeneracy. The Chern 
number is always an integer, and for a simple double de- 
generacy, it equals ±1. In other words, near a DP, B„ 
has the form of a monopole field with flux equal to ±27r. 
Hence, to find the DPs, we numerically evaluate B„ for 
an initial H, and follow the lines of B„ in the direction 
of increasing strength until we hit a monopole. Since the 
number of DPs where levels n and n + 1 are degenerate is 
topologically fixed and known, all the DPs can be found 
by taking sufficiently many initial values of H. The DPs 
for successive pairs of levels occur in layers, with H z es- 
sentially constant in a layer. With h = g/isH/2fciJ, 
the first three layers are at h z = (exactly), 0.0427, and 
0.0854. It should be noted, however, that in a given layer, 



(m, m') 


Layer 


(h x ,h y ,h z ) 


(-10, 10) 


1 


(0.0404, 0, 0) 
(0.1207,0,0) 
(0.2011,0,0) 
(0.2565,0,0) 
(0.3139, ±0.0642,0) 
(0.3875, ±0.1334,0) 
(0.4656, ±0.2043,0) 


(-10,9) 


2 


(0,0,0.0426) 
(0.0787, 0, 0.0426) 
(0.1568,0,0.0427) 
(0.2276,0,0.0427) 
(0.2893, ±0.0447, 0.0427) 
(0.3594, ±0.1109, 0.0427) 
(0.4340, ±0.1788, 0.0427) 


(-10,8) 


3 


(0.0385,0,0.0853) 
(0.1150,0,0.0853) 
(0.1892,0,0.0853) 
(0.2650, ±0.0268, 0.0854) 
(0.3326, ±0.0900, 0.0854) 
(0.4038, ±0.1548, 0.0854) 


(-9,9) 


1 


(0.0384,0,0) 
(0.1147,0,0) 
(0.1889,0,0) 
(0.2651, ±0.0261,0) 
(0.3325, ±0.0891,0) 
(0.4034, ±0.1535,0) 


f-8 8) 


1 


(0.0365, 0, 0) 
(0.1089,0,0) 
(0.1805,0,0) 
(0.2339,0,0) 
(0.2831, ±0.0511,0) 
(0.3468, ±0.1092,0) 



one can have DP's corresponding to tunneling between 
levels with different pairs of Zeeman quantum numbers. 
For example, layer 1 contains points corresponding to 
tunneling between m = -10 <H> 10, -9 «-> 9, — 8 <H> 8, 
etc. Likewise, layer 2 contains points corresponding to 
tunneling between m = -10f>9, —9 <->■ 8, —8 f-> 7, etc. 
In Table [H we show the DP's for (-10 O 10) (-10 O 9), 
(-10 O 8), (-9 O 9), and (-8 O 8) tunneling. For the 
first three of these pairs of levels, the projections of the 
DP's onto the xy plane are shown in Fig. Q] We note 
that it is just these three pairs of levels for which tunnel 
splittings were reported in Ref. Q. Hoewever, except for 
the ground pair of levels (layer 1), not even all the DP's 
on the x axis are found there. 

In the rest of this note, we discuss the form of the Berry 
curvature near a DP in more detail. For simplicity, we 



3 



will divide H by kiJ 2 . Since V H % and E n are both 
divided by this factor, it follows from Eq. ([2]) that B„ is 
unchanged. With this preliminary remark, let us suppose 
that E n = E n _i = E nn _i at h = h„ . n _i, and denote 



r = h — h,. 



(4) 



Further, let us make a particular choice of the two de- 
generate states at r = 0, and denote them by |o) and \b), 
with (b\a) = 0. (Any orthogonal linear combination of 
| a) and \b) would also work.) It suffices to truncate the 
Hamiltonian to this two dimensional subspace since the 
sum in Eq. ([2]) is dominated by degenerate states. Hence, 
at r = 0, we have 



rl — E n n -i 



1 0' 
1, 



(5) 



For small enough r, we can take \a) and \b) to be un- 
changed, so 



^_ I Jaa Jab 
J \ J&q Jbb , 



(6) 



where J aa = (a|J|a) etc. Next, let us define J aa + J;,;, = 
Jt, J aa - J bb = Ju, 3 ab = J(v + w)/2, where t, u, v, 
and w are real vectors. In terms of these vectors, we have 



H = — 



t ■ r / u • r (v + iw) • r 
J \(v-jw)-r u ■ r 



, (?) 



where we have ignored the constant E n . n -i. Similarly 
ignoring the overall shift — t • r/J, the eigenvalues of this 
matrix are ±e(r), with 

e(r) = [(u • r) 2 + (v • r) 2 + (w • r) 2 ] 1 / 2 . (8) 

To write the eigenvectors compactly, we define 

cos6»(r) = u-r/e(r), (9) 
sin6»(r)e^ (r) = (v + iw) ■ r/e(r). (10) 



The eigenvectors are then 

, 1+) 



sin £fe 



-up 



sin 



(11) 



Further abbreviating c = cos ^8 and s = sin ^6, and 
g = (+|VhK|— }, we have 

g = -2csu - s 2 e~ %v (v + iw) + c 2 e lv (v - iw). (12) 

It then follows that 

g x g* = 2i[cos6>(v x w) + sin#cos!y2(w x u) 

+ sin0siny>(u x v)], (13) 



so that for the level labeled 
-1 



2e 3 (r) 

-1 
2e 3 (r) 



[(u • r)(v x w) + (v • r)(w x u) 

+ (w • r)(u x v)] 
[u • (v x w)] r. (14) 



This is clearly of monopole form with appropriately 
scaled and sheared axes. It is not difficult to show that 
B_ = -B+ and that V • B+ = 0. 

To find the Chcrn number, we must evaluate the inte- 
gral 



Q ± = --L 4 B. .,/ ;i s 



(15) 



for a suitable surface S. Let us take S to be the surface of 
the parallepiped with vertices at (±u* , ±v* , ±w* ) , where 
u*, v*, and w* are the reciprocal vectors 



v x w w x u 
V = w 

n ' n ' 

fi = u • (v x w). 



(16) 
(17) 



Then u ■ u* = 1, u • v* = u • w* = 0, etc. Let us consider 
the integral over the face of the parallepiped that has 
edges along v* and w*, and thus has a normal parallel 
to +u (the others have normals along — u, ±v, and ±w). 
On this face, we may parametrize r as 



r = u* + av* + /3w*, 



(18) 



where — l<a<l, — 1</3<1. Then, since the area of 
this face is 4|v* x w*| = 4|u|/0, 



d s= ^da d(3. 



(19) 



Finally, u • r = 1, v • r = a, and w • r = /3. Thus, the 
contribution from this face to Q± is given by 



47T 



1 r l 



da d/3 1 
(l + a 2 + /3 2 ) 3 / 2 = 6' 



(20) 



(The integral is elementary, and can be performed by 
standard trigonometric substitutions.) The contributions 
from the other faces are identical, so 



Q-i 



±1. 



(21) 



Thus the Chern number associated with the lower en- 
ergy level is always —1 irrespective of the details, and this 
is why the sum of the numbers for a given energy level 
cannot be altered by varying the anisotropy parameters. 

The reader cannot have failed to notice the regularity 
of the DP pattern for the present problem. This is in 
sharp contrast to the generic situation where no pattern 
is expected apart from a general scaling of the density 
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with energy |3| . The reason for this regularity is unclear 
to us. 

This work is supported by the NSF via grant no. PHY- 
0854896. 
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